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^ ; ABSTRACT 

O I We present new results of a spectroscopic survey of circumstellar Hi in the 

I direction of evolved stars made with the Nangay Radiotelescope. The Hi line 

I at 21 cm has been detected in the circumstellar shells of a variety of evolved 

. ^ [ stars : AGB stars, oxygen-rich and carbon-rich. Semi- Regular and Miras, and 

^ I Planetary Nebulae. The emissions are generally spatially resolved, i.e. larger 

I than 4', indicating shell sizes of the order of 1 pc which opens the possibility 

to trace the history of mass loss over the past ~ lO'* — 10^ years. The line- 
profiles are sometimes composite. The individual components have generally a 
quasi- Gaussian shape; in particular they seldom show the double-horn profile 
that would be expected from the spatially resolved optically thin emission of a 
uniformly expanding shell. This probably implies that the expansion velocity 
decreases outwards in the external shells (0.1-1 pc) of these evolved stars. 

The H I line-profiles do not necessarily match those of the CO rotational lines. 
Furthermore, the centroid velocities do not always agree with those measured in 
the CO lines and/or the stellar radial velocities. The Hi emissions may also 
be shifted in position with respect to the central stars. Without excluding the 
possibility of asymmetric mass ejection, we suggest that these two effects could 
also be related to a non-isotropic interaction with the local interstellar medium. 
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Hi was detected in emission towards several sources (p Per, a Her, Lyr, 
U CMi) that otherwise have not been detected in any radio hnes. Conversely 
it was not detected in the two oxygen-rich stars with substantial mass-loss rate, 
NML Tau and WX Psc, possibly because these sources arc young with hydrogen 
in molecular form, and/or because the temperature of the circumstellar Hi gas 
is very low (< 5 K). 

Subject headings: stars: AGB and post-AGB - (stars:) circumstellar matter - 
stars: late-type - stars: mass-loss - (ISM:) planetary nebulae: general - radio 
lines: stars 



1. INTRODUCTION 

Hydrogen is the most abundant element in the atmospheres of Asymptotic Giant Branch 
(AGB) stars and in their outflows. It may be found in the form of atomic hydrogen. Hi, 
or molecular hydrogen, H2, besides relatively minor (in abundance) species. The fraction 
of hydrogen in each of these two forms is an important factor controlling the atmospheric 
structures of red giants and also the winds which develop from their atmospheres (e.g. Lamers 
& Cassinelh 1999). 

It is therefore of interest to study the circumstellar environments of AGB stars in the 
low excitation hydrogen lines. Unfortunately molecular hydrogen is presently difficult to 
detect because its lowest rotational lines lie in a region of the spectrum which cannot be 
accessed from the ground, and because a high spectral resolution is needed. 

On the other hand atomic hydrogen has a line at 21 cm in an easily accessible, and 
well protected, region of the radio spectrum. However the observations are difficult be- 
cause the intensity is weak and above all because there is a contamination by interstellar 
hydrogen along the same lines of sight. For a long time these difficulties have limited the 
Hi observations of AGB sources to o Get (Bowers & Knapp 1988, BK1988). Nevertheless 
we recently readdressed this question with the upgraded Nangay Radiotelescope and could 
detect this line in emission towards several AGB sources. We have already presented our 
results on four of them: RS Gnc (Gerard & Le Bertre 2003, Paper I), EP Aqr and Y GVn 
(Le Bertre k Gerard 2004, Paper II) and X Her (Gardan et al. 2006, Paper III). We de- 
tected these four AGB sources with a good signal-to-noise ratio and showed that the line 
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profiles are composite. The individual components have Gaussian-hke profiles, although in 
some cases a rectangular profile could be fitted as well. The H I emissions are spread over a 
velocity range confined to the total CO velocity extent. However, we noted that sometimes 
the central velocities in CO and H I do not coincide exactly. Furthermore, in Paper 1, we 
argued that atomic hydrogen is present already in the stellar atmosphere of the central star 
(RS Cnc) which might even be devoid of molecular hydrogen. In Paper II we could resolve 
angularly two circumstellar shells, and show that H I may be present at large distances from 
the central stars, 0.5 pc or more for EP Aqr. The Hi data display a large scale asymmetry 
suggesting that asymmetric structures may show up early in the AGB phase (EP Aqr & X 
Her), well before the planetary nebula phase. In the case of the source associated to Y CVn, 
the H I data even seem to probe a large region extending out to where the stellar outflow is 
interacting with the Interstellar Medium (ISM). 

It thus appears that the H I emission line profiles, when observed with a good signal-to- 
noise ratio and a fair removal of the ISM contamination, may reveal important clues on the 
physics of stellar outflows and their kinematics. Also, as atomic hydrogen has a long lifetime 
in the ISM, it can be used to trace circumstellar envelopes out to large radial distances. In 
the present paper, we give results on 22 new sources for which we have now good quality 
spectra, or good upper limits. The main emphasis of our programme is put on AGB sources 
and we have tried to select a variety of representative specimens in terms of variability types, 
chemical composition, mass loss rates, etc. In addition, we have also included two planetary 
nebulae in order to sample more completely the post-main-sequence evolution of stars of low 
and intermediate masses. Two of our sources have already been observed in Hi with the 
Very Large Array (VLA), o Ceti (BK1988) and NGC7293 (Rodriguez et al. 2002, R2002). 
This allows a useful comparison between single-dish and interferometric data. 



2. OBSERVATIONS 

2.1. Observing Procedure 

The observations have been performed with the Nangay Radiotelescope (NRT). The 
NRT is a meridian instrument with a rectangular effective aperture of 160 m x 30 m. The 
HPBW at 21 cm is 4' in the East- West direction and 22' in the North-South direction. It 
has been upgraded (van Driel et al. 1996) and reopened to observations in 2001. The point 
source efficiency is 1.4 K Jy~^ at 21 cm and the total system noise ~ 35 K. The autocorrelator 
has 8192 channels split in 4 banks of 2048, recording the linear (PA = and 90°) and circular 
polarisations. For our programme we consider only the total intensity (Stokes parameter I). 
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Two modes of observation have been used. In the frequency-switch (f-switch) mode, the 
data are acquired on source with a spectral resolution of 0.08 kms~^ and a bandwidth of 160 
kms~^. This mode is very useful to check the level of interstellar Hi emission in the direction 
of the source. A first indication of this level can be obtained by consulting the existing 
galactic Hi surveys (e.g. Hartmann & Burton 1997), but their sensitivity limits are too high 
to readily select directions in which it might be possible to detect circumstellar Hi with 
a priori minimum background confusion. In the course of our programme, we characterize 
quahtatively the confusion by "weak" , "mean" and "high" . These indications are reported 
in Table 2 (col. 2, see Sect. 2.2). They are based on Hartmann & Burton's results that are 
obtained with a single dish of HPBW ~ 0.5°, as well as on our own f-switch observations. 
It should be noted that the impact of confusion depends on the telescope beam size and on 
the observing strategy, in particular single-dish versus interferometric mode. Except ionnally 
Hi from the source can be seen directly on the f-switch spectra (e.g. Y CVn, Paper II). 

For the circumstellar H I observations proper we use the position-switch mode, in which 
the telescope field of view (f.o.v.) is alternated between the on-source position and 2 off- 
source positions, placed symmetrically in the East- West direction. We refer to Paper II for 
a more detailed description of our observing procedure. It has proven to be most effective 
in separating the circumstellar emission from the interstellar one. A reason for this success 
comes from the fan beam shape (4' x 22') of the telescope that allows a good East- West 
discrimination of the circumstellar envelope against the galactic H I background. It should be 
noted that this procedure works perfectly when the intensity of the background varies linearly 
in right ascension across the region centered on the source. However, if the background 
presents a quadratic variation, or variations of higher even orders, spurious spectral features 
may appear. Such features generally grow rapidly with increasing separation between the 2 
off-positions, which allows to identify them. 

In general we selected a bandwidth of 80 kms~^ with a spectral resolution of 0.04 kms~^. 

For some sources where a large expansion velocity is expected from CO line data or where 
galactic confusion affected part of the spectrum, we have used a bandwidth of 160 kms~^. 
The data are then smoothed to give a final spectral resolution of 0.32 kms~^. Except when 
mentioned otherwise, this resolution is used by default in the figures (Sect. 2.2). 

Finally, the source sizes are often comparable or larger than the NRT beam (4' in the 
East- West direction). The intensity increases with beam throw and then reaches a maximum 
for ± n beams when both off-positions (East and West) lie entirely outside the source. We 
then reconstruct the total intensity (stot) by summing up the source fluxes in the different 
off-positions (s^) and in the central position {sq): 
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n-l 

Stot = XI «fe (1) 
k=-n+l 

This is done by combining the position-switch spectra obtained at ± k beams (CkEW) in 
the following expression: 

n— 1 

Stot = (2n - l)xCnEW - 2 ^ W (2) 

fc=i 

n being such that the intensity in 2 successive position-switch spectra is the same, i.e. 
Cn+iEW — CnEW. The radius of the source can then be estimated to (n — 1) x 4' ± 2'. 
Equ. (2) assumes that the background varies linearly across the sky. On the other hand it 
makes no assumption on the brightness distribution of the source. For an unresolved source 
(0 < 4'), Equ. (2) reduces to: 

Stot = CiEW = C2EW = C3EW, etc. (3) 

In the figures (next Section), the position-switch spectra {CiEW, C2EW, etc.) are shown in 
the upper panels, and their averages and the final spectra (stot) are displayed at the bottom 
(continuous and dashed curves, respectively). 

The average spectrum has a physical meaning only for an unresolved source (Equ. 3). 
However, the resulting signal-to-noise ratio is in general higher than for the total spectrum, 
so that spectral details are sometimes better seen on the average profile (see, e.g. below, Y 
UMa or R Peg). The difference between the two sets of spectra is the result of the source 
resolution by the 4' beam. 

2.2. Sample of Sources 

Table 1 lists the sources discussed in the present paper. We give the common name used 
in this work (col. 1), the IRAS identifier (col. 2), and some basic data. Except if mentioned 
otherwise, the spectral type (col. 3), the variabihty type (col. 4) and the period (col. 5) 
are taken from SIMBAD, directly or through the VizieR Service. The distances (col. 6) 
are in general deduced from the Hipparcos parallaxes. Next (col. 7), we give an estimate 
of the stellar effective temperature, Teff . This datum is important because hydrogen in the 
stellar atmosphere and in the inner circumstellar envelope is expected to be mostly atomic 
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for Tefr> 2 500 K and molecular for Tefr< 2 500 K (Glassgold & Huggins 1993, GH1983). 
Finally we give the LSR radial velocity (Vigr, col. 8), the expansion velocity (Vexp, col. 9), 
and the mass loss rate, estimated in general from CO rotational line data. These data are 
useful because both the central velocity and the velocity extent of the H I lines arc expected 
to be close to those given by the CO lines. In general CO data from different observers are 
consistent within ±2 kms~^. When no CO data are available (5^ Lyr, p Per, Her, U 
CMi, NGC 6369), we adopt equivalent parameters from optical observations. Unfortunately, 
optical data seem less rehable than CO data (±5 kms~^), perhaps in part due to the stellar 
variabihty. Except for U CMi, circumstellar envelopes have been detected around all sources 
of our sample, either in the CO radio lines, or in optical lines. 

In Table 2, we give the parameters characterizing the Hi emissions. In general these 
emissions can be fairly well fitted with a Gaussian. We derive the LSR central velocity, Vcent. 
(col. 4), a full width at half maximum, FWHM (col. 5), and a peak intensity, Fpeak (col. 6), 
obtained from an unconstrained positive Gaussian fit to the integrated flux (stot; Equ. 2). 
In the cases of composite profile ( Y UMa, R Peg) , two such Gaussians are also fitted to the 
average profile in the lower panel of the corresponding figure. When part of the spectrum is 
affected by interstellar confusion it has been excluded from the range used for the Gaussian 
fit. In such cases, the part of the spectrum which is affected is always that closest to 
kms~^ LSR. The Flux (col. 7) is estimated from the integral of the Stot Gaussian fit and 
therefore only takes into account the extension of the source in right ascension. Finally the 
hydrogen mass, Mhi (col. 8), is derived from the Flux and for the distance, d (Table 1). 
We assume negligible absorption of the circumstellar H I emission by foreground galactic gas 
and self- absorption, and thus may slightly underestimate the H I mass. 



2.2.1. AGB Sources 

Stars on the AGB present a large variety of properties in terms of effective temperature, 
luminosity, chemistry, variability, etc., probably reflecting sources of different ages, initial 
masses, etc. The observed mass loss rate ranges from 10~^ to 10""^ Moyr~^. In fact AGB 
stars seem to undergo mass loss at a rate which may vary on timescalcs ranging from a few 
years to perhaps 10^ years. We have tried to select oxygen-rich and carbon-rich sources with 
different variability types and different mass loss rate estimates. Our sample covers also a 
wide range of stellar temperatures. 
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Table 1. Observed sources in right ascension order and basic data. 



Source 


IRAS name 


Spectral 


Variability 


Period 


d 


Tcff 




Vcxp 


M 






type 


type 


(days) 


(pc) 


(K) 


(km s^-"-) 


(kms^-"-) 


(M0yr-i) 


WX Psc 


01037+1219 


M8-M10: 


Mira 


660 


650 


<2500 


+9 


23 


2 10-5 


o Get 


02168-0312 


M7IIIc+Bcp 


Mira 


332 


128 


~2700 


+46 


7 


2.5 lO-'^ 


p Per 


03019+3838 


M4n 


SRb 


40 


100 


3576 


+25 


<8.5 


1 10"* 


NML Tau 


03507+1115 


M6e-M10e 


Mira 


462 


245 


2000-3000 


+35 


19 


3 10-s 


S CMi 


07299+0825 


M6e-M8e 


Mira 


332 


365 


2500-3100 


+52 


3 


2 10-'^ 


U CMi 


07386+0829 


M4e 


Mira 


414 


191 


3400 


+42 






RV Hya 


08372-0924 


M5II 


SRc 


116 


318 


3420 


-43 


5 


3 10-^ 


U Hya 


10350-1307 


C5II 


SRb 


450 


162 


2965 


-31 


7 


9 10"** 


Y UMa 


12380+5607 


M7II-III 


SRb 


168 


313 


3100 


+19 


5 


2 10-7 


RY Dra 


12544+6615 


G4,5J 


SRb 


200 


488 


2810 


-4 


13 


1.5 lO-'^ 


RT Vir 


13001+0527 


M8ni 


SRb 


155 


138 


3034 


+17 


9 


1 lo-''' 


W Hya 


13462-2807 


M7.5c-M9ep 


SRa 


361 


115 


2500 


+41 


8 


2 10-'^ 


Her 


17123+1426 


M5Ib-n 


SRc 




117 


3285 


-17 


10 


1 10-7 


NGC 6369 


17262-2343 


PN 






1550 


> 10000 


-90 


41 




52 Lyr 


18527+3650 


M4n 


? 


60-120 


275 


3460 


-7 


<10 


5 10-8 


Z Cyg 


20000+4954 


M5e-M9e 


Mira 


264 


490 


<3300 


-148 


4 


4 10-8 


NGC 7293 


22267-2102 


PN 






200 


> 10000 


-23 


32 




RPeg 


23041+1016 


M6e-M9e 


Mira 


378 


350 


2300-2900 


+24 


9 


5 10-7 


AFGL 3068 


23166+1655 


C 


Mira 


696 


1140 


1800-2000 


-31 


14 


1 lO-'' 


AFGL 3099 


23257+1038 


C 


Mira 


484 


1500 


1800-2000 


+47 


10 


8 10-s 


TX Psc 


23438+0312 


C5II 


Lb 




233 


3115 


+13 


10 


2 10-7 


R Gas 


23558+5106 


M7HIe 


Mira 


430 


107 


<3000 


+25 


12 


5 10-7 



-8- 



Oxygen-rich AGB Sources 

We discuss first 5^ Lyr, whose variability type is not defined, and then the 6 oxygen-rich 
Semi- Regular variables of our sample. 

(5^ Lyr (HR 7139) is a relatively warm AGB star (M411), which is barely variable (A V ~ 
0.2 mag., Bakos & Tremko 1991). We adopt the effective temperature determined by Sudol 
et al. (2002) from near-infrared (J and H) interferometry, 3 460 ± 160 K. From the profiles 
of photospheric lines affected by circumstellar absorption, Sanner (1976) finds evidence of an 
outflow. This source has not been detected in any radio fine, but it has an IRAS extended 
counterpart (Young et al. 1993a). On the other hand, the ISO (3-40 iim) spectrum (Heras 
et al. 2002) shows no evidence of dust. In Table 1, we adopt the parameters of its wind from 
Sanner (1976). 

The Hi spectra are presented in Fig. 1. In this figure, and in the following ones, the 
horizontal bar represents the velocity extent (Visr± Vexp) expected from the wind parameters 
given in Table 1 (full line in the case of a CO radio-hne determination, dotted line in the 
case of an optical determination). The contamination is serious especially at V > — 5 kms~^. 

However, each spectrum clearly shows an Hi emission that is increasing with offset up to 
=b 16'. This indicates that the source has a diameter of ~ 24'. The central velocity of the 
profile and the FWHM agree with those obtained by Sanner (1976). The intensity has been 
integrated over the 7 NRT beams that contribute to the signal. At a distance of 275 pc, 
it translates to ~ 0.075 M© in atomic hydrogen. As the source has a size of ~ 24' in the 
East- West direction and as the NRT beam is 22 ' in the North-South direction, we may have 
missed part of the Hi emission (see Sect. 3.2 for a corrected Hi mass). 
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-50 50 

LSR velocity (km/s) 

Fig. 1. — Top panel: 5^ Lyr spectra obtained in the position-switch mode. The off-positions, 
in this figure and ensuing ones, are indicated in bold characters above each spectrum. For 
clarity, the individual spectra have been displayed with vertical offsets of 0.05 or 0.1 Jy. 
Bottom panel: average of these 5 spectra and space integrated intensity scaled by a factor 
1/2. The horizontal dashed bar indicates the velocity range expected for circumstellar matter 
(Sanner 1976). 
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Fig. 2. — Top panel: p Per spectra obtained in the position-switch mode. Bottom panel: 
average of the 5 spectra and space integrated intensity scaled by a factor 1/5. The horizontal 
dashed bar indicates the velocity range expected for circumstellar matter (Sanner 1976). 
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p Per (HR 921) has not been detected in CO, nor in any other radio hne. Sanner (1976) 

finds evidence of an outflow with parameters (Visr, Vexp, M) given in Table 1. Mauron & 
Guilain (1995) detected circumstellar Nal at 5" from the central star (i.e. ~ 500 AU), 
confirming that the star is losing matter. However, the infrared spectrum obtained by 
ISO shows no evidence of dust around p Per (Heras et al. 2002). We adopt the effective 
temperature from Dumm & Schild (1998). 

An Hi emission is barely detected in the spectrum obtained with off-positions at ±4', 
whereas it is clearly detected at ±16', and at ±24' with basically the same intensity (Fig. 2). 
This indicates that the H I emission is resolved and has a size corresponding to about 7 NRT 
beams in the East- West direction, i.e. ~ 24'. Galactic confusion clearly affects the Hi profile 
near 20 kms~^ at ± 16 and 24'. The blue wing of the Hi profile is cut down, which probably 
results in a displacement of the H I velocity centroid. 
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Fig. 3. — Top panel: a Her spectra obtained in the position-switch mode. Bottom paneh 
average of these 5 spectra and space integrated intensity scaled by a factor 1/5. The hor- 
izontal dashed bar indicates the velocity range expected for circumstellar matter (Deutsch 
1956). 
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Her (HR 6406) is known to undergo mass loss since Deutsch (1956) showed that 
circumstellar Unes are present in absorption in the spectrum of its companion, at an angular 
distance of 5". We adopt his wind parameters. Like p Per, it was not detected in any radio 
line (e.g. Heske 1990). Circumstellar Kl and Nal were also detected in emission by Mauron 
& Caux (1992). We adopt the effective temperature determined by Perrin et al. (2004a) 
from 2 fim interferometric data (Tefr = 3 285 ± 89 K). 

The Hi emission is increasing with offset (Fig. 3), which indicates that the source has 
a diameter of the order of 16'. The central velocity of the profile and its width agree with 
the parameters obtained by Deutsch (1956). The interstellar confusion is high for V > —5 
kms~^, and seems even to affect the red wing of the profile obtained at ±16'. An emission 
component stands out at —10 kms^^. As it can be seen on the individual position-switch 
spectra at ±4', ±8', and ±12' it is probably real. 

RV Hya has been detected in CO by Winters et al. (2003). We adopt their parameters 
for the wind characteristics. With the temperature scale of normal giants (Ridgway et al. 
1980) and an M5 spectral type, the effective temperature should be around 3400K. 

The Hi line is marginally detected on the average (Fig. 4). It does not seem to be 
angularly resolved. 
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LSR velocity {kr\/s) 



Fig. 4. — Average spectrum of RV Hya obtained in the position-switch mode. The spectral 
resolution is 0.64 kms~^. 
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Y UMa is associated to an IRAS extended source (Young et al. 1993a). It has been 
observed in CO (2-1) by Knapp et al. (1998) and their wind characteristics are given in 
Table 1. These estimates agree with those obtained by Olofsson et al. (2002). From the 
temperature scale of normal giants of type later than M6 (Perrin et al. 1998), we adopt an 
effective temperature of about 3 100 K. 

The Hi emission is spatially extended (0 ~ 8', Fig. 5). Wc note that the central 
velocity of the Hi feature is shifted by ~ 2 kms~^ with respect to the CO emission and is 
narrower (Knapp et al. 1998). The profile is probably composite with a narrow component 
ovcrimposcd on a broad faint one. A detailed examination of the separate position-switch 
data shows that the Hi emission is asymmetric and slightly shifted West. 
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Fig. 5. — Top panel: Y UMa spectra obtained in the position-switch mode. Bottom panel: 
average of these 5 spectra and space integrated intensity. Gaussian fits to the average are 
shown to illustrate the decomposition of the profile in 2 components. 
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RT Vir is also associated to an IRAS extended source (Young et al. 1993a). We adopt 
the wind characteristics from Knapp ct al. (1998). However, Olofsson et al. (2002) derive 
a larger mass loss rate estimate (~ 5xl0~^ MqYt'^, at 170 pc). It has been detected in 
the thermal (v=0, J=2-l) SiO line by Gonzalez Delgado et al. (2003) at a slightly different 
velocity (Visr~ 18.6 kms^^). This source was also detected in the OH maser lines at 1612, 
1 665 and 1 667 MHz by Etoka et al. (2003). The effective temperature is taken from Dumm 
& Schild (1998). 

The Hi source is extended, (p ~ 24' (Fig. 6). A detailed examination of the position- 
switch data shows that the H I emission is asymmetric and shifted East by ~ 4'. 
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W Hya is sometimes classified as a Mira (e.g. Young 1995). It is associated to an IRAS 
extended source (Young et al. 1993a). We adopt the wind characteristics from Knapp et al. 
(1998). These parameters agree with those obtained by Young (1995) and Olofsson et al. 
(2002). The H2O line at 557 GHz was detected by the Odin satellite and gives consistent 
parameters for the wind (Justtanont et al. 2005). The circumstellar shell was also detected 
in Kl fluorescence by Guilain & Mauron (1996). Like RT Vir, it has been detected in the 
OH maser lines at 1612, 1665 and 1667 MHz (Etoka et al. 2003). An effective temperature 
of 2 500 ± 190 K was obtained by Haniff et al. (1995) from optical interferometry. The 
stellar disk was also resolved in the radio continuum by Reid & Menten (1997) who derived 
a (radio) brightness temperature of 1500 ±570K. 

The H I source is clearly resolved with a maximum flux reached only at ±16' so that we 
estimate its size to ~ 24' (Fig. 7). The spectra obtained at ±16' and ±24' show a narrow 
dip at +20kms~^ (FWHM ~ 3 kms~^), reveahng the presence of a small interstellar cloud 
at ~ 20' East of W Hya. 
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50 
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Fig. 6. — Top panel: RT Vir spectra obtained in the position-switch mode. Bottom paneh 
average of these 5 spectra and space integrated intensity scaled by a factor 1/2. 
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We now discuss the 8 oxygen-rich Miras of our sample, in order of increasing spectral 
type, which corresponds approximately to increasing period and increasing observed mass 
loss rate. The spectral types of Miras are known to vary with phase, and accordingly the 
stellar effective temperatures, which makes them particularly uncertain. 

U CMi is probably one of the bluest Miras. As it has not been detected in any radio 
lines, we adopt the radial velocity determined optically (General Catalogue of Stellar Radial 
Velocities, GCRV). Adopting the effective temperature scale for Mira variables of Dyck et 
al. (1974), Teff should be around 3400 K. 

As there is no estimate available for the expansion velocity, the stellar radial velocity 
is only marked with a vertical bar (Fig. 8). The Hi line is easily detected with a triangular 
profile centered at a velocity (40.9 kms^^) very close to the stellar one (42 kms~^). It is 
more intense for ±8' than for ±4' and we estimate its diameter to ~ 8'. As for Y UMa, we 
observe that the Hi emission is not centered on the star position, and slightly shifted West. 
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Fig. 7. — Top panel: W Hya spectra obtained in the position-switch mode. Bottom paneh 
average of these 6 spectra and space integrated intensity scaled by a factor 1/5. 
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Z Cyg is remarkable for its high radial velocity suggesting it might belong to population 

II. It has been detected as an OH maser by Sivagnanam et al. (1989) and in the CO (3-2) 
line by Young (1995) from whom we adopt the outflow parameters and the distance. The 
effective temperature should be around 3 300K or below (Dyck et al. 1974). 

An Hi hne is suspected at ~ -146 kms~^ (Fig. 9) but we only set an upper hmit to its 

flux. 
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Fig. 8. — Top panel: U CMi spectra obtained in the position-switch mode. Bottom paneh 
average of these 4 spectra and space integrated intensity. 
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S CMi has been observed in CO (1-0) and (2-1) by Winters et al. (2003). We adopt 
their wind parameters and distance based on the period-luminosity relationship for Miras 
(Feast 1996). The effective temperature should stay most of the time above 2 500 K (Dyck 
et al. 1974). 

The Hi hne is easily detected (Fig. 10) and is shghtly more intense for ±8' than for ±4'. 
As there is no clear difference between ±8' and ±12' the diameter of the source should be 
~ 8'. 
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Fig. 9. — Average spectrum of Z Cyg obtained in the position-switch mode. The spectral 
resolution is 0.64 kms~^. 
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o Get (Mira, HR 681) has been observed extensively. Its infrared spectrum shows a 
lO-fim sihcate feature strongly in emission (e.g. Speck et al. 2000). The CO profiles are 
complex, possibly composite (Knapp et al. 1998; Winters et al. 2003). Circumstellar Kl 
and Nal were also detected by Mauron & Caux (1992). CO (2—1) interferometric and Kl 
data show that the envelope is strongly asymmetric, perhaps bipolar, (Josselin et al. 2000). 
We adopt a mass loss rate of 2.5x10"^ M© yr~^ from the detailed modelling of CO rotational 
lines by Ryde & Schoier (2001). The central star has a hot companion, at approximately 
0.5" (~ 65 AU), surrounded by a disk accreting a part, «i 1 %, of the matter in the outflow 
(Reimers & Cassatella 1985). 

On the basis of a large set of infrared interferometric data, Weiner (2004) finds that a 
2 200 K H2O shell surrounds a 2 700 K star (at a phase close to maximum). This estimate of 
the stellar temperature agrees with the Dyck et al. (1974) scale. On the other hand, Perrin et 
al. (2004b), using near- infrared interferometric data, derive an effective temperature in the 
range 3200-3 600 K. These results illustrate that Mira atmospheres have complex structures, 
as predicted by hydrodynamical models (e.g. Woitke et al. 1999). A prediction of the atomic 
hydrogen abundance in the stellar atmosphere and inner outflow would certainly require a 
detailed hydrodynamical and chemical modeling. 

o Cet was already detected in Hi with the VLA (BK1988). These early observations 
suggest that approximately one third of the hydrogen is atomic when it leaves the star. The 
source is spatially resolved (radius ~ 160") and slightly offset to the North- West. However 
this large-scale asymmetry does not seem to be related to the bipolar flow (0 ~ 20") observed 
in CO by Josselin et al. (2000). 

The signal that we obtained is slightly larger for ±8' than for ±4' and very stable for 
larger offsets (Fig. 11). We estimate that the source is just resolved spatially with the NRT 
(0HI ~ 8 ± 4'), consistent with the result of BK1988 and the extension (0 ~ 4.4') found 
by IRAS at 60 fim (Young et al. 1993a). The average line-profile is of high quality and 
resembles the triangular shape of the CO line-profiles (e.g. Knapp et al. 1998; Winters et 
al. 2003), with a small asymmetry shifting the peak (46.3 kms~^) by ~ 1 kms~^ to the 
red of the centroid velocity (45.7 kms~^). Such a close agreement, in shape and in position, 
is exceptional in our sample. But the most important fact maybe that this high-quality 
H I line-profile clearly does not show the rectangular shape expected from an optically thin 
envelope with constant expansion velocity (Le Bertre & Gerard, Paper II). This point is 
further discussed in Sect. 3. 
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Fig. 10. — Top panel: S CMi spectra obtained in the position-switch mode. Bottom paneh 
average of these 3 spectra and space integrated intensity. 
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The integrated intensity is ~ 0.56 Jyxkms""^ which (for a distance of 128 pc) translates 
into Mhi~ 2.2x10"^ Mq. For comparison, BK1988 find 0.42 Jyxkms"^, which is in good 
agreement with our estimate. 

For a size of 8' (0.30 pc) and an average expansion velocity of ~4 kms~^, the crossing 
time is ~ 38x10^ years, and the average mass loss rate in Hi over this period is ~ 0.6x10"^ 

M0yr~^. Assuming atomic hydrogen represents 3/4 of the total mass, wc find (M)av ~ 
0.8x10^^ MQyr^^ which is one third of the present mass loss rate estimated from CO-lines 
modelling (Ryde & Schoicr 2001). This result is in line with the interpretation of BK1988 
that 1/3 of the hydrogen is leaving the star in atomic form. 

R Peg has been observed in CO by Young (1995) and Winters et al. (2003). Wc adopt 
their wind characteristics, van Belle et al. (1996) obtained interferometric data at 2 ^m 
and derived effective temperatures at phases close to maximum (2 333 ± 100 K), and close 
to minimum (2 881 ± 153 K). 

As for Y UMa the H I profile seems composite with a narrow component overimposed on 
a broader one (Fig. 12), and is reminiscent of the Y CVn Hi profile (Paper 11). Furthermore, 
the narrow component (1) is spatially resolved with a size of ~16'. However, the difference 
in the central velocities of the 2 components is noteworthy. The central velocity of the 
spectrally broad component (2) seems to agree with the CO one found by Winters et al. 
(2003), +24 kms"^, and with those determined from the OH masers: +23.9 kms~^ for 1667 
MHz and +24.6 kms"^ for 1665 MHz (Sivagnanam et al. 1989). 

The examination of the separate position-switch data shows evidence of an offset West 
with respect to the central star. 
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Fig. 11. — Top panel: o Get spectra obtained in the position-switch mode. Bottom panel: 
average of these 5 spectra and space integrated intensity. 
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R Cas (HR 9066) has been observed in CO by Young (1995) and Knapp et al. (1998). 
We adopt their wind characteristics. Neri et al. (1998) have obtained maps in the CO (1-0) 
and (2-1) hnes and find the source to be extended, ~ 17" and 21", respectively. It is 
also associated to a 60 /xm IRAS extended source of 8.6' (Young et al. 1993a). From 2 fim 
interferometric data van Belle et al. (1996) determine an efi^ective temperature of 2 954 ± 
174 K, at a phase close to maximum. 

The Hi source is barely detected at ±4' (Fig. 13), but clearly at ±12'. The diameter is 
thus ~16'. 

NML Tau (IK Tau) is a large amplitude Mira whose spectral type exhibits important 
variations with phase (from M6 to MIO). It has been detected in CO (1—0) by Knapp 
& Morris (1985) and Nyman et al. (1992), and mapped in the CO (2-1) line by Neri et 
al. (1998, (p ~ 17"). This source presents a slightly self-absorbed silicate feature at 10 /jm 
(Speck et al. 2000). Adopting the scale for Mira variables of Dyck et al. (1974), the effective 
temperature should vary in the range 2 000-3 000 K. We adopt the period from the infrared 
monitoring of Le Bertre (1993), the distance from the period-luminosity relation for 0-rich 
Miras (Le Bertre & Winters 1998) and the mass loss rate derived from the modelling of the 
infrared energy distribution (Le Sidaner & Le Bertre 1996). It was not detected in Hi by 
Zuckerman et al. (1980). 

Our Hi data suffer from interstellar confusion over the velocity range —40, -|-40 kms~^. 
Fig. 14 shows no convincing signal that could be attributed to NML Tau. Over the remaining 
positive velocity range defined by the CO fine- profile (-1-40, -|-54 kms~^) we find an upper 
limit to the flux density level of 20 mJy. or 0.8 Jyxkms"^ over the full CO velocity extent 
of 38 kms^^. At 245 pc, this limit translates to Mhi < 0.01 Mq in the volume covered by 
the "on" position of the NRT beam. 

This limit applies to an Hi emission line. However, some circumstellar matter could 
be at a low temperature such that atomic hydrogen would appear in absorption against the 
background (see Sect. 4.3). An upper limit to the background level is given by the sum of 
the interstellar Hi emission around +35 kms~^ (~ 0.5 K) and of the continuum emission at 
1420 MHz (~ 3.5 K, Reich & Reich 1986), therefore - 4 K. 
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Fig. 12. — Top panel: R Peg spectra obtained in the position-switch mode. Bottom panel: 
average of these 4 spectra and space integrated intensity scaled by a factor 1/2. Gaussian 
fits to the average are shown to illustrate the decomposition of the feature in 2 components. 
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Fig. 13. — Top panel: R Cas spectra obtained in the position-switch mode. Bottom paneh 
average of these 5 spectra and space integrated intensity scaled by a factor 1/2. 
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WX Psc (IRC +10011) is the most extreme oxygen- rich AGB source in our sample and 
is a representant of the type II OH/IR class, although at a high galactic latitude (b = —50°). 
It has been detected in CO (1—0) by Knapp & Morris (1985). The CO emission has been 
found to be extended by Neri et al. (1998, ~ 20"). The spectral type is probably strongly 
variable, as for NML Tau. The IRAS Low Resolution Spectrum shows a clear self-absorbed 
10 /xm silicate feature (Volk & Cohen 1989). The effective temperature should stay most of 
the time under 2 500K (Dyck et al. 1974). As for NML Tau, we adopt the period from 
Le Bertre (1993), the distance from Le Bertre & Winters (1998) and the mass loss rate from 
Le Sidaner & Le Bertre (1996). Zuckerman et al. (1980) also attempted to detect Hi at 21 
cm, but without success. 

Our data suffer from interstellar confusion up to -|-10 kms~^ (Fig. 15). As for NML Tau 
there is no convincing signal that can be attributed to WX Psc. In the region > 10 kms~^, 
we can set an upper limit of 5 mJy to the Hi emission. Following the same procedure as 
for NML Tau, we find an upper limit of 0.23 Jyxkms"^ to the Hi flux within the central 
beam. At 650 pc, this hmit translates to Mhi < 0.02 M© in the volume probed by the "on" 
position of the NRT beam. 

The same caveat as for NML Tau apphes if this upper limit is interpreted as an absorp- 
tion line, the background in that direction being at a level of 4.5 K. 



-34- 




-50 50 100 

LSR velocity (km/s) 



Fig. 14. — Average of all the spectra obtained in the position-switch mode on NML Tau. The 
spectral resolution is 0.64 kms"^. The horizontal bar indicates the velocity range expected 
from the CO emission. 
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Fig. 15. — Average of all the spectra obtained in the position-switch mode on WX Psc. The 
spectral resolution is 0.64 kms~^. The horizontal bar indicates the velocity range expected 
from the CO emission. 
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Carbon-rich AGB Sources 



We first discuss one Irregular and two Semi-Regular variables, then two high mass loss 
Miras, sometimes referred to as extreme carbon stars (ECS). As in the oxygen-rich cases, 
the effective temperatures of carbon Miras are variable and very uncertain, particularly since 
these two are infrared sources. 
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Table 2. Gaussian fits to the H l profiles, integrated intensities, and atomic hydrogen 

masses. 



Source 


confusion 


integration 


Figure 


Vcent. 


FWHM 




Flux 


Mhi 






time (hours) 




(kms-1) 


(kms-i) 


(mjy) 


(Jyxkms~l) 


(Mo) 


WX Psc 


higii/weak 


30 


15 






< 5 


< 0.23^^ 


< 0.02'' 


Get 


weak 


34 


11 


+45.1 


7.7 


73 


0.56 


0.0022 


p Per 


medium 


50 


2 


+30.3 


10.6 


193 


2.05 


0.0048 


NML Tau 


high/ weak 


31 


14 






< 20 


< 0.76*^ 


< O.Ol'' 


S CMi 


weak 


27 


10 


+51.0 


5.2 


44 


0.23 


0.0072 


U CMi 


mean 


20 


8 


+41.3 


8.1 


119 


0.96 


0.0083 


RV Hya 


weak 


35 


4 


-43.8 


8.5 


6 


0.05 


0.0012 


U Hya 


weak 


42 


17 


-30.0 


4.9 


102 


0.50 


0.0031 


YUMa 


weak 


28 


5 


+16.5 


3.7 


98 


0.36 


0.0084 


RY Dra 


weak 


25 


18 


-4.5 


8.5 


407 


3.46 


0.194 


RT Vir 


mean/ weak 


24 


6 


+13.2 


9.6 


192 


1.84 


0.0083 


W Hya 


high / mean 


43 


7 


+41.3 


21.8 


549 


12.0 


0.037 


Her 


mean/high 


46 


3 


-14.4 


11.1 


206 


2.29 


0.0074 


NGC 6369 


mean 


22 


21 


-91.8 


40.5 


31 


1.26 


0.71 


52 Lyr 


mean/high 


21 


1 


-10.3 


5.7 


737 


4.20 


0.075 


ZCyg 


weak 


27 


9 






<6 


< 0.05*^ 


< 0.003=- 


NGC 7293 


weak/mean 


34 


22 


-28.7 


34.8 


530 


18.44 


0.174 


R Peg 


weak 


30 


12 


+23.5 


5.3 


166 


0.88 


0.025 


AFGL 3068 


high 


23 


20 


-26.8 


30.9 


80 


2.47 


0.76 


AFGL 3099 


weak 


49 


19 






< 2 


< 0.04" 


< 0.02'" 


TX Psc 


high 


24 


16 


+9.8 


4.2 


1250 


5.25 


0.067 


R Cas 


mean 


33 


13 


+24.5 


5.5 


122 


0.67 


0.0018 



''in the central beam ( "on-source" ) 
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TX Psc (HR 9004) has a CO envelope with an irregular structure (Heske et al. 1989). 

The CO rotational line profiles may be composite with a narrow peak and a broad component 
(sec also Heske 1990). Wc adopt the wind parameters from Loup et al. (1993) and the 
effective temperature (3115 ± 130 K) from Bcrgcat et al. (2001). The source has been 
found to be extended at 60 /im by IRAS (Young et al. 1993a). 

The Hi signal grows rapidly with increasing offset (Fig. 16). We estimate the Hi diam- 
eter at 24'. 
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Fig. 16. — Top panel: TX Psc spectra obtained in the position-switch mode. Bottom paneh 
average of these 5 spectra and space integrated intensity scaled by a factor 1/5. 
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U Hya (HR 4163) has been detected in CO (2-1) by Knapp et al. (1998) and we adopt 
their wind parameters. The IRAS data show a detached dust shell (0 ~ 3.5') produced by 
a previous episode of much higher mass loss rate, ~ 5x10"^ Moyr'^ (Waters et al. 1994). 
We take the effective temperature from Bergeat et al. (2001). 

The Hi line is faint (Fig. 17), but well detected on the average spectrum, and might be 
composite. Also the emisson is spatially resolved {(f) ~ 32') . Although the flux is weak, this 
case might be similar to Y CVn. 

RY Dra belongs to the rare J-type class, hke the star Y CVn which was discussed in 
Paper II. It has been detected in the ^^CO and ^^CO rotational lines by Jura et al. (1988), 
the latter with about half the intensity of the former. The ISO data show a slight evidence 
of a detached shell that Izumiura & Hashimoto (1999) ascribe to a previous episode of mass 
loss at a rate of ~ 10~^ MqJT''^. Furthermore, the IRAS data show a very extended source 
(0 ~ 40') at 60 and 100 nm (Young et al. 1993a). We adopt the effective temperature from 
Bergeat et al. (2001). 

The Hi source is weakly detected at small offsets (Fig. 18) but reaches more than 50 
mJy at ±24'. The diameter is the largest in our sample, ~ 40'. 
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Fig. 17. — Top panel: U Hya spectra obtained in the position-switch mode. Bottom paneh 
average of these 6 spectra and space integrated intensity scaled by a factor 1/2. 
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The 2 following sources were discovered in the AFGL infrared survey. Their carbon-rich 
nature was established from near-infrared spectrophotometry by the presence of the 3.1 //m 
carbon-star absorption feature (Gehrz et al. 1978; Jones et al. 1978). 
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Fig. 18. — Top panel: RY Dra spectra obtained in the position-switch mode. Bottom paneh 
average of these 7 spectra and space integrated intensity scaled by a factor 1/5. 
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AFGL 3099 (IZ Peg) was detected in CO (1-0) by Knapp & Morris (1985). The CO 

(2—1) emission has been mapped by Neri et ah (1998) and shown to be asymmetric (14" x 7"). 
We adopt the period from Le Bertre (1992), the distance from the period-luminosity rela- 
tionship for carbon Miras (Groenewegen & Whitelock 1996) and the mass loss rate from the 
modelling of the infrared energy distribution (Le Bertre 1997). There is no stellar spectral 
classification; however, as the 3.1 /xm feature is extremely deep (Gehrz et al. 1978), which is 
an indication of a very low temperature, we adopt the effective temperature for the coolest 
carbon stars in Bergeat et al. (2001). 

Wc could not detect the source for any beam-throw (up to ± 6 beams, i.e. ±24'). In 
Fig. 19 we present the average of all our H I spectra. This is the deepest integration in our 
survey. In the 30-80 kms~^ range the H I interstellar background is low and perfectly removed 
by our procedure, so that we can safely set an upper limit of 2xl0~^ Jy to any circumstellar 
H I emission feature in the central position. Adopting the CO expansion velocity, we find an 
upper limit of 0.04 Jyxkms"^ to the Hi emission, which, at a distance of 1.5 kpc, translates 
to Mhi < 2xl0~^ Mq in the volume probed by the on-position. 

AFGL 3068 (LL Peg) has such a dense wind that the SiC dust feature at 11 yum is seen 
in absorption (Jones et al. 1978). It has been detected in CO (1—0) by Knapp & Morris 
(1985) and Nyman et al. (1992). The CO (2—1) emission has been mapped by Neri et al. 
(1998) who suggest the presence of a large (0 ~ 17") detached bipolar shell. Mauron & 
Huggins (2006) have detected galactic light scattered by circumstellar dust out to a distance 
of ~ 40" from the central star with an amazing spiral structure. We adopt the period and 
the distance from Le Bertre (1992, 1997). A mass loss rate of 5x10^^ MqJt^^ was derived 
from a radiative transfer model of the infrared energy distributions by Le Bertre et al. 
(1995) whereas a consistent time-dependent hydrodynamical model of this source obtained 
by Winters et al. (1997) gives ~ 1x10"^ Mgyr-^ As for AFGL 3099, we adopt the effective 
temperature for the coolest carbon stars (Bergeat et al. 2001). 

An emission feature centered at —31 kms~^ is barely present at ±4', but clearly detected 
at ±8', ±12' and ±16' (Fig. 20). The intensity is increasing from ±4' to ±8', but not 
significantly from ±8' to ±16'. We conclude that H I is detected in emission from AFGL 3068 
and that this emission is extended with a size of at least ~ 8' in diameter, or 2.6 pc at 1.14 
kpc. 
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Fig. 19. — Average of all the spectra obtained in the position-switch mode on AFGL 3099. 
The spectral resolution is 1.28 kms~^. The horizontal bar indicates the velocity range 
expected from the CO emission. 



-46- 



2.2.2. PN Sources 

After they have expelled most of their stellar envelope red giants are thought to evolve 
toward the Planetary Nebula (PN) stage. The mass loss decreases and the hot central core 
progressively emerges and starts to ionize its surrounding. For some time, a neutral remnant 
from the AGB circumstellar shell can still be present around an ionized nebula. Recent 
developments may imply that this scenario applies only to the giants with more massive 
progenitors, and/or to those that are in a close binary system and experience a common 
envelope interaction. 

"Circumnebular" neutral hydrogen has been detected in absorption against the con- 
tinuum emitted by the central ionized nebula in 5 PNs by Taylor et al. (1990). Atomic 
hydrogen has also been detected in emission from IC 418 (Taylor et al. 1989) and from NGC 
7293 (R2002). 

NGC 6369 (Little Ghost Nebula) is a PN of optical diameter ~ 40". The central 
star (HD 158269) is an hydrogen-deficient [WC] Wolf-Rayet star. A detailed photoionization 
model of the nebula by Monteiro et al. (2004) shows that it has a hourglass structure. 
They derive a nebular mass of 1.8 M©. We adopt the distance derived from this modelling 
and the expansion velocity (Vcxp= 41.5 kms~^) determined from the [O III] line width by 
Meatheringham et al. (1988). We note that the [O III] line should probe the ionized central 
part of the nebula, and that it probably overestimates the expansion velocity of the external 
medium. Indeed, Garay et al. (1989) have obtained an WIQa linewidth of only 42.8 kms~^. 
This source has been seen to be extended by IRAS at 25 //m (0 ~ 1.1'); this extension could 
be related to the un- identified 21-//m emission feature reported by Hony et al. (2001). 

The source is barely detected at ±4' (Fig. 21, but clearly at ±8'. We estimate its 
diameter at 8'. 
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Fig. 20. — Top panel: AFGL3068 spectra obtained in the position-switch mode. Bottom 
paneh average of these 4 spectra and space integrated intensity. 
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NGC 7293 (Helix Nebula) is a nearby PN (~200pc) with a large extent on the sky 
(0 ~1000"). Young et al. (1999) have mapped the nebula in the CO (2-1) line. The 
CO linewidth is about 63 kms^^. The emission is localized in condensations tracing an 
expanding equatorial ring that is seen almost pole-on and that surrounds the bright ionized 
nebula (Huggins ct al. 1999). Speck et al. (2002) have obtained H2 (v=l— 0) line map which 
shows that this emission also coincides with the bright nebula. Furthermore they present an 
Ho; image which shows a weak and more extended emission, up to ~ 1100" from the central 
star. This weak emission coincides with dust emission detected by IRAS at 60 and 100 
lim, and by ISO at 160 fj,m. Neutral carbon (C I hne at 492 GHz) was detected in several 
positions on the bright nebula by Young et al. (1997). 

NGC 7293 was observed in Hi with the VLA (R2002). The emission is fragmented and 
follows the CO (2—1) ring, except in the South-East quadrant where it is missing. With a 
diameter of ~ 16', it surrounds the bright optical nebula; in particular it is absent in the 
central part (0 ~500"). R2002 estimate the Hi mass to ~ 0.07 Mq. However, from the 
continuum data at 21 cm, they estimate that they may have missed ~ 50% (0.67/1.4) of 
the flux. 
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Fig. 21. — Top panel: NGC 6369 spectra obtained in the position-switch mode. The spectral 
resolution is 0.64 kms~^. Bottom panel: average of these 3 spectra and space integrated 
intensity. The horizontal dashed bar indicates the velocity range expected from the [O III] 
line- width (Meatheringham et al. 1998). 
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The interstellar contamination seems large at V > —20 kms~^. The source is clearly 
detected in emission for V < —20 kms~^ (Fig- 22) and throws larger than 8'. We detect 
nearly the same intensity at ±16' and ±24', and thus evaluate the diameter at ~ 24'. This 
is much larger than estimated by R2002, but consistent with the finding of Ho; and dust 
emissions up to 18' from the central star by Speck et al. (2002). 

A comparison of our Hi spectrum with the spatially integrated spectrum of R2002 
shows that the narrow feature at —15 kms^^ probably belongs to the source and is not due 
to confusion, i.e. the contamination may affect our data only at V > —15 kms^^. Also, the 
intensity of this feature levels off at ± 12', an indication that it is not due to the galactic 
background. Nevertheless we have derived Hi line parameters (Table 2) by using only the 
part of the spectrum outside the range (—20, ±20 kms~^). We find a central velocity ~ —29 
kms~^, instead of —24 kms~^ (Young et al. 1999), probably because we excluded the —15 
kms~^ feature from the fit. Also we observe in our average spectrum a narrow component 
at —50 kms~^ that is not seen in the R2002's spectrum. However, a CO feature is detected 
at this velocity by Young et al. (1999). The latter is concentrated in the South-East part of 
the Helix ring where R2002 detect no H I emission. 

For the main component centered at —29 kms~^, the examination of the separate 
position-switch data does not show evidence of a spatial offset with respect to the cen- 
tral star. On the other hand the emission at —15 kms^^ is offset to the West. R2002 have 
mapped this emission with the VLA (their figure 7). The corresponding channel map shows 
that the feature is located in the North- West and that it could be due to an interaction with 
the local ISM. 

For the main component, we measure an integrated Hi intensity of 18.4 Jyxkms^^, 
which is a factor 2.3 larger than the R2002's estimate. This is in agreement with their finding 
that, with the VLA, they may have lost ~ 50% of the flux. The H I emission missed by R2002 
is probably diffuse and extends beyond the classical torus of diameter 16'. The intensity 
that we missed in the narrow feature at —15 kms~^ is only 0.5 Jyxkms"^. Considering 
the extent that we find and the North-South size of our beam (22'), it is fikely that we also 
missed H I fiux, at the level of ~ 30 %, to the North and to the South of the central source, 
bringing the total fiux to ~ 28 Jyxkms"^. A full map similar to those that we did in Paper 
H for EP Aqr and Y CVn would be useful. 

The Hi fiux from NGC 7293 is indeed so large that the source can be seen on the 
Hartmann & Burton (1997) maps obtained at velocities close to —40 kms~^ for which the 
galactic confusion is moderate. 



-51 - 




-100 -50 50 

LSR velocity (km/s) 



Fig. 22. — Top panel: NGC 7293 spectra obtained in the position-switch mode. The spectral 
resolution is 0.64 kms~^. Bottom panel: average of these 5 spectra, and space integrated 
intensity scaled by a factor 1/5. 
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3. INTERPRETATION 
3.1. Line Profiles 

One of the most striking results emerging from our survey is that we find H I line-profiles 
with a quasi-Gaussian shape, whereas we were expecting rectangular profiles, for unresolved 
sources, or double-horn ones, for resolved sources. Sometimes the profile seems triangular 
with a peak close to the centroid velocity (e.g. RT Vir, Fig. 6, or o Get, Fig. 11). 

In Paper II, we indeed performed numerical simulations of the H I emission produced by 
a spherical circumstellar envelope. The emission is assumed to remain optically thin, and its 
brightness temperature to be proportional to the H I column density. Therefore there is no 
shell self- absorption and the velocities are radial. In these conditions for an unresolved source 
with a uniform wind velocity and a constant production rate in atomic hydrogen, the line- 
profile is rectangular, is centered on the stellar radial velocity, and has a width equal to twice 
the expansion velocity. For a resolved source, the line profile shows two horns with peaks 
separated by twice the expansion velocity. These simulations were performed assuming a 
uniform response within the telescope beam. However, introducing a more realistic response 
did not significantly change the results (see e.g. figure 7 in Paper III) . 

We seldom observe double-horn profiles. Only a Her may show such a shape. We 
also note that NGG7293 shows several components that could be interpreted as horns. 
Likewise, the profiles of AFGL 3068 and NGG 6369 could be non-Gaussian. Previously, in a 
few peculiar cases, we detected a component of the H I emission that might be fitted by a 
rectangle (RS Gnc, Paper I; EP Aqr, Paper II; X Her, Paper HI). In these cases, the width 
is relatively narrow, ~ 3-4 kms~^. 

To get centrally peaked profiles we had to assume that the H I velocity is varying with 
distance from the central star. For resolved sources the velocity has to decrease to zero 
relative to the stellar radial velocity. A large quantity of hydrogen has to accumulate close 
to zero velocity, otherwise a deficit appears in the center of the profile. For unresolved 
sources, the velocity could be assumed to decrease as well as to increase with distance. 
However, it would be unlikely to get increasing velocities only for unresolved sources, and 
therefore we can conclude that, on scales probed by the Hi emission, the expansion velocity 
always decreases with distance from the central star. This slowing-down of the circumstellar 
wind could be due to its interaction with external matter, either remains of previous episodes 
of mass loss and/or surrounding ISM (Young et al. 1993b; Paper II). 
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Fig. 23. — Half-FWHM of the Hi line versus expansion velocity determined from CO rota- 
tional lines. 
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The slowing-down is also suggested by the comparison of the expansion velocity derived 
from the CO line-profiles to the half-FWHM of the H I line (Fig. 23). Except for AFGL 3068 
and W Hya, the former is always larger than the latter. We note that AFGL 3068 is at a 
large distance from the galactic plane (~ 740 pc), perhaps in a region devoid of interstellar 
matter. 

In general a gaiissian profile can result from thermal and turbulent broadening, as well 
as from systematic motions of the gas. Presently in our modelling, we are not taking into 
account thermal and turbulent effects. On the one hand their inclusion in the model may 
force to reduce the H I outflow velocity because they will provide a supplementary source 
of line broadening. On the other hand the velocity gradient may have to be reduced and 
outflow velocities close to (relative to the central star) may not be needed (Libert, Gerard 
& Le Bertre, in preparation). 



3.2. Spatial Distributions 
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Although the spatial resolution of the NRT is low, most of our sources are so large that 
they are resolved in the East- West direction. In Table 3 we give estimates of the diameters 
of the Hi sources (col. 3) as explained in Sect. 2.1. For sources that are not resolved with 
the NRT, we adopt a diameter of 2'. The four sources that were not detected in H I emission 
(WX Psc, NML Tau, Z Cyg and AFGL 3099) have not been reported in this Table. We 
then correct (col. 4) the measured Hi mass assuming that the Hi emission is circularly 
symmetric, i.e. it has the same size in declination and in right ascension. The sizes of the 
Hi shells (col. 5) are given for the distances adopted in Table 1. A characteristic timescale 
(col. 6) can be derived from this size and the FWHM of the H I emission line. Finally from 
the estimated mass of hydrogen in the circumstellar shell, and assuming that 75 % of the 
mass is in hydrogen, we derive an average mass loss rate over this timescale (col. 7). 
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Table 3. Estimated characteristics: sizes, timescales, average mass loss rates. 
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Fig. 24. — Hi versus IRAS 60 //m diameter from Young et al. (1993a). 
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Young et al. (1993a) measured the angular extent of red giants and young PNs in the 
IRAS 60 and 100 /im survey data. For comparison, we give also the diameter estimated by 
these authors from the 60 ^m IRAS data (col. 8). In general we find sizes in Hi that are 
comparable or larger (Fig. 24). This should not be overemphasized as Hi diameters are only 
rough estimates. 

o Get was barely resolved with the NRT, although it is one of the closest sources in 
our sample. We find an estimate of the size of the Hi source (in the East-West direction) 
which is consistent with the size of the dust shell estimated from the IRAS data at 60 fim 
(R ~ 2.2', or 0.08 pc. Young ct al. 1993a). Wc note that BK1988 detected Hi emission up 
to ~ 160" (2.7') from the central star. At the other extreme, AFGL3068, which is at ~ 1.2 
kpc, was angularly resolved with the NRT, ~ 8', or > 2 pc. The H I circumstellar shells 
of Mir as may therefore have very different physical sizes, suggesting very different durations 
for the past mass loss. 

We also note that in the case of AFGL 3068 the H l shell seems much larger than the shell 
observed in the far-infrared that was not resolved by IRAS (Young et al. 1993a) or by ISO 
(Izumiura & Hashimoto 1999). This may be an effect of the dust emission dependence on 
temperature that may impede the detection of matter far from the central star or indicate 
a relatively recent onset of copious mass-loss rate that is not reflected in the HI spectra. 
Interestingly, our H I data imply a mass-loss rate about 10 times lower than that given by 
the detailed modelling of Winters et al. (1997). 

ISO observations at 60 and 90 yum of WX Psc show a point source with a FWHM ~ 1' 
(Hashimoto & Iziimiura 2000). For a distance of 650 pc, the circumstellar shell diameter 
should be at most 0.2 pc. This is interpreted as a recent evolution with a sudden event 
(< 4000 years) of extremely high mass-loss rate. However, Hashimoto & Izumiura do not 
exclude the possibility of a faint extended emission that would be the trace of a past mild 
mass-loss phase. A recent onset of the mass loss in sources like WX Psc or NML Tau (that 
was not resolved by IRAS, Young 1993a) might explain their non-detection in Hi. Indeed 
in young compact sources with central star effective temperature lower than 2500 K, most 
of the circumstellar matter may be molecular (see Sect. 4.4). 

We have also examined our position-switch data by using separately the two off-positions. 
The profiles are severely affected by the spatial variations of the background and the resulting 
spectra suffer much more from galactic H I confusion, which makes them in general difficult 
to exploit. In some cases (weak/mean confusion, large signal), e.g. RT Vir, U GMi, we can 
detect the circumstellar emission and observe that it is not symmetric with respect to the 
central position (i.e. the star position). Nevertheless, this should be taken with caution as 
the increasing confusion towards kms~^ velocity may affect the baselines of the separate 
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East and West profiles. However, the same effect was noted previously for EP Aqr (Paper 
II) and X Her (Paper III) for which we also observed shifts in the North-South direction. 
Finally wc note that in the o Ceti map obtained with the VLA by BK1988 the Hi emission 
appears offset to the North- West. 

Prom IRAS data at 60 //m Young et al. (1993b) find that among nearby stars Semi- 
Regular variables are more likely to be resolved than Miras. Our statistics is small but 
our results (at least for oxygen-rich sources) tend to agree with their conclusion that Semi- 
Regular variables have been losing mass for a longer duration than Miras. On the other 
hand they do not find evidence that the shapes of the shells arc distorted by interaction 
with the ISM, which we do (see also Sect. 4.2). Finally, AGB carbon-rich sources appear 
more extended than oxygen-rich ones. The difference is still more striking when comparing 
the total Hi masses (col. 3 in Table 3). This quantity ranges from 0.001 to 0.1 Mq for 
oxygen-rich sources, and from 0.005 to 0.7 M© for carbon-rich sources. 

The larger extension of carbon-rich sources as compared to oxygen-rich ones is not 
so surprising as carbon-rich AGB stars should be in a phase of evolution that follows the 
oxygen-rich stage. And the longer the star loses mass, the larger the mass of its circumstellar 
shell. 



3.3. H I versus CO 

We have detected H l emission in sources with and without emission from the CO rota- 
tional hues. For instance a}- Her is detected with a peak intensity ~ 50 mJy at ±12', but was 
not detected in CO by Heske (1990) although the presence of circumstellar matter close to 
the central star was demonstrated by Deutsch (1956). We note that this situation appears 
preferentially for warm giants with Teff> 3 000 K. It suggests that CO is underabundant in 
these sources for which the outflows might be mostly atomic. If this suggestion is correct 
then these sources could be detectable in O I and C I lines, for instance at 492 GHz. Knapp 
et al. (2000) did not detect Cl at 492 GHz from Her, but this could be an effect of a 
cyclic variation of the mass loss rate, the region that they probed being only 15" in diameter. 
We also note that our estimate for the average mass loss rate is 2.5 times larger than that 
of Deutsch (1956). 

We sometimes easily detect Hi from Semi- Regular variables with low mass loss rates 
whereas we may not detect it from Miras with large mass loss rates. This may be an effect of 
the dependence of the H2/H1 equilibrium in the stellar atmosphere on temperature. Indeed 
the stellar effective temperature of the former are larger than 2 500 K, and smaller for the 
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latter, and GH1983 show that hydrogen should be mostly atomic for Teff> 2 500 K, and 
molecular for Teff< 2 500 K. Atomic hydrogen should then be present only in the external 
parts of the circumstellar shells of Miras with large mass loss rates where H2 is photo- 
dissociated by the Interstellar Radiation Field. This point is further discussed in Sect. 4.3 
and 4.4. CO and H I can trace outflows only where these species are present. 

If we now compare CO and H I profiles of sources that show both emissions, wc find that, 
with the notable exception of o Get, they often differ. CO profiles are generally parabolic 
or rectangular whereas H I profiles are Gaussian or triangular. The H I emission is found 
within the range encompassed by the CO emission. The Hi FWHM is generally smaller 
than the expansion velocity determined from CO data (Fig. 23). Sometimes its centroid is 
shifted by 1-3 kms"^ with respect to the CO centroid. For instance, in Y UMa, we find 
Hi emission from 14 to 20 kms~^ with a centroid at 17 kms~^, whereas the CO emission 
spans from 14 to 24 kms~^ and is centered at 19 kms~^ (Knapp et al. 1998). Although the 
kinematics of the outflows should be the same in CO and H I, the regions that are probed 
by these two tracers are different. Firstly, the beams are different, typically 10-50" for CO 
observations, and 4'x22' for our Hi observations. Secondly, CO and atomic hydrogen may 
be present in different zones of the circumstellar shells. For instance CO is expected to be 
photodissociated by the interstellar radiation field at distances ~ 10^^-10^^ cm depending 
on the mass loss rate (Mamon et al. 1990), whereas Hi should in general be protected from 
photo-dissociation by the surrounding ISM. Therefore CO and Hi trace different parts of 
the circumstellar shells. Finally the CO hne-profiles and intensities are expected to depend 
on excitation and optical depth effects much more than for H I. 
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When we compare the mass loss rates estimated from the Hi hne profiles to those 
obtained from CO rotational lines, we do not find any clear trend (Fig. 25). This probably 
indicates that the mass loss rate varies by an order of magnitude over timescales ~ 10^ years. 
Nevertheless, we note that Miras tend to lie below the one-to-one line, and Semi-Regulars 
above it. 

It thus appears that the H I line at 21 cm and the CO rotational lines arc complementary 
tracers of the circumstellar outfiows and that they are equally needed to describe these media. 

4. DISCUSSION 

4.1. Hi Confusion 

In general our data are affected by the large scale structure of the galactic H I emission. 
Thus the final limits of detection are not set by the sensitivity of the observations but by 
the effectiveness in removing the galactic emission. When the position-switch procedure has 
been effective we have quahfied the confusion as weak (Table 2) . For "medium" and "high" 
confusion part of the Hi spectrum can be affected. It indicates that matter in a direction 
close to the target and at velocities in the expected range has been detected. In general 
this matter, that can be at any distance along the line of sight, should not be related to the 
stellar source. 

However, as the central stars may have been undergoing mass loss for a long duration 
and/or during a previous stage of evolution (e.g. the RGB), in some cases it might be that 
this matter (or part of it) is physically related to the source. For instance, the interaction 
of the outflow with its surrounding leads to the formation of a shell made of compressed 
circumstellar material and swept-up circumstellar /interstellar matter (Sect. 3.1). The matter 
in this shell should be at velocity close to the central star velocity. If part of the material lost 
by the central star is mixed in position and in velocity with the surrounding ISM, we may 
never be able to get an exhaustive balance of the mass lost by the central stars. For instance, 
in their Hi study of IC 418, Taylor et al. (1989) noted that part of the detected hydrogen, 
blueshifted with respect to the PN system velocity, cannot be identified unambiguously as 
circumstellar or interstellar matter. 
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4.2. Interaction with the ISM 

For several sources we observe a shift between the central velocities of the H I and CO 
lines (~ 1-3 kms~^). In some cases we also observe an offset of the Hi source with respect 
to the stellar position (~ 1-4'). Such a kind of offset has not been observed in CO. Since 
the Hi emission profiles show evidence of a slowing- down of the stellar outfiows by the 
surrounding ISM (Paper II and HI, Sect. 3.1), and since the central stars are moving with 
respect to the ISM^, we may suspect that these velocity shifts and spatial offsets are due to a 
non-isotropic interaction between the outfiows and the ISM. One may also consider a gradient 
in the ISM density that would induce a non-isotropic slowing-down of the circumstellar 
outfiow. 

In order to corroborate this hypothesis we have used the non-spherical model of H I emis- 
sion developed by Gardan et al. (2006, Paper HI). As an example we consider a source lo- 
cated at 140 pc with an outflow corresponding to a mass loss in atomic hydrogen of 5xl0~^ 
M0yr^^, and that is initially spherical and isotropic. The velocity is taken to decrease 
linearly from an inner boundary located at 0.1' (or 4x10"^ pc) from the central star. To 
simulate a non-isotropic interaction, we adopt a velocity law as : 

V = 10.0 - (1.43 + QAlcosa) x (r - 0.1) (4) 

with V expressed in km s~^, r in and a being a polar angle. The slowing down of the outfiow 
is maximum for a = 0° and minimum for a — 180°. The geometry of the circumstellar shell 
is thus progressively becoming egg-shaped. For this test, we adopt a typical duration of 
50x10^ years, such that the outer boundary is located at 5' (0.2 pc) for a — 0°, with V = 
1 kms-\ and 7.2' (~ 0.3 pc) for a = 180°, with V ~ 2.8 kms-^ We adopted for the NRT 
beam a "sine" response (see Paper HI). 

The results of such a simulation are given in Fig. 26. The polar axis of the source {a — 
0°) is first placed in the plane of the sky pointing East (Fig. 26, left). The Hi profiles are 
symmetric and centered on the stellar velocity, but the emission feature is more pronounced 
East than West. If the polar axis is pointing 90° away from the observer (Fig. 26, right), the 
H I profile becomes asymmetric with a red-shifted peak. In that case the emission remains 
centered on the star position. Finally, if the polar axis is oriented 45° away from the observer 
and at a position angle of 90° (Fig. 26, center), the Hi profile is asymmetric and the emission 
is stronger East of the central star. 



^although this cannot be strictly demonstrated, we note that our sources are close to the Sun and have 
a non-zero Vigr 
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This simulation shows that a slowing-down of the stellar outflow more pronounced in a 
given direction could produce either a shift in velocity of the H I emission line, or an offset 
in position of this emission, or a combination of both. This effect may add to an intrinsic 
asymmetry due to the mass loss from the central star, as we showed for X Her (Paper III). 

4.3. Cold Atomic Hydrogen 

We have detected Hi in emission towards AFGL3068. However, we have previously 
reported the detection of Hi in absorption towards a similar source, IRC +10216, (Lc Bertre 
& Gerard 2001), indicating that the gas in the outflow cools down to a temperature below 
4 K. Such a low temperature was predicted for high mass loss rate AGB stars as a result 
of adiabatic expansion (Jura ct al. 1988; Sahai 1990). At first sight it is puzzling that 
we do not observe the same effect in AFGL 3068 of mass loss rate about 5 times larger. 
However the distances to these sources differ by a factor ~ 8, so that the NRT beam scales 
as 0.16pcx0.86pc for IRC+10216 as compared to 1.3pcx7.3pc for AFGL3068. Therefore, 
in the case of AFGL 3068 we have explored a very external region where grain photo-electric 
heating might become important. For instance, in their models of IRC+10216, Crosas & 
Mentcn (1997) and Groenewegen et al. (1998) predict that the temperature should increase 
at r > 0.1 pc due to this effect (although their minima are still higher than our observational 
upper limit). 

The Hi detection of IRC+10216 in absorption and of AFGL 3068 in emission might 
then be explained as a combined effect of beam size and distance. If this is correct, the 
Hi mass that we have derived for AFGL 3068 may be underestimated due to absorption 
against the background. 

The gas temperature in the circumstcllar envelopes depends strongly on the mass-loss 
rate (e.g. Jura ct al. 1988). It is expected to be larger in the sources with low mass-loss rate 
(~ 10^^ M0yr^^) than in those with high mass-loss rate (~ 10~^ MQyr"^). Low mass-loss 
rate sources, and in particular Semi-Regular Variables, are therefore less liable to contain 
very cold atomic hydrogen in their envelopes. On the other hand, this effect may concern 
Miras with substantial mass loss rate, hke WX Psc or NML Tau. 

4.4. Molecular Hydrogen 

Hydrogen in stars with low atmospheric temperature (T^ < 2500 K) should be molecular 
(GH1983). The expected lifetime of molecular hydrogen in the ISM is ~ 1000 years for a 



-64- 



standard Interstellar Radiation Field (Le Petit et al. 2002). This translates to r ~ 10~^ or 
10~^ pc for an average expansion velocity of 1 or 10 kms~^, respectively. However, molecular 
hydrogen is self-shielded (Morris & Jura 1983) and the real extent of an H2 circumstcllar 
envelope depends on the mass loss rate and the velocity field. For instance GH1983 estimate 
that, for a mass loss rate of 4xl0~^ M© yr^^ and a uniform expansion velocity of 16 kms~^, 
the radius of the molecular envelope would be ~ 0.2 pc. 

It is therefore important to check the presence of atomic hydrogen at large distance 
from the central star for those with low effective temperature. As sources with low effective 
temperature are also those with presently high mass- loss rate (e.g. WX Psc, NML Tau), 
atomic hydrogen may be detected in these sources, either in absorption or in emission, only 
at large distances from the central stars, and, of course, only if mass loss has been operating 
for a long enough time. 

5. CONCLUSION 

We have undertaken a search for the H I line at 21 cm in the direction of AGB stars and 
related sources. We have detected circumstcllar H I in a variety of sources sampling the post 
main-sequence evolution of low and intermediate mass stars. The emissions are extended 
indicating shell sizes from 0.2 pc to more than 2 pc. The H I hne can therefore be used to 
trace the mass loss history for long periods, up to several 10^ years. However, we may have 
missed atomic hydrogen at low temperature (< 10 K) and hydrogen locked into H2. 

For some sources (p Per, a Her, P Lyr, U CMi) our Hi detection is the first radio 
detection. 

The Hi line-profiles indicate that the outflows are slowing down, probably due to the 
interaction with surrounding matter. Furthermore, we observe offsets in velocity and spatial 
asymmetries that indicate a distortion of the outflows. These effects were already noted for 
EP Aqr (Paper II) and for X Her (Paper HI), and seem frequent. They could be connected 
to the motion of the sources with respect to the local ISM, or to its inhomogeneity, as well 
as to the mass ejection process itself. 

We are grateful to the referee for his/her detailed review that helped us improve our 
paper. The Nancay Radio Observatory is the Unite scientifique dc Nancay of the Obser- 
vatoire de Paris, associated as Unite de Service et de Recherche (USR) No. B704 to the 
French Centre National dc la Recherche Scientifique (CNRS). The Nancay Observatory also 
gratefully acknowledges the financial support of the Conseil Regional de la Region Centre in 
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Fig. 25. — Average mass loss rate estimated from H I versus mass loss rate from the literature 
(Table 1). Filled circles: mass loss rates derived from CO rotational lines; empty circles: 
mass loss rates derived from optical data (Deutsch 1956; Sanner 1976). 
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Fig. 26. — Simulation of the Hi emission from an egg-shaped circumstellar shell (see 
Sect. 4.2). Left : PA = 90°, inchnation angle = 0°; centre : PA = 90°, inchnation an- 
gle = -45°; right : PA = 90°, inchnation angle = -90°. 



